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Abstract Computational modelling is helpful for elucidat¬ 
ing the cellular mechanisms driving biological morphogene¬ 
sis. Previous simulation studies of blood vessel growth based 
on the Cellular Potts model (CPM) proposed that elongated, 
adhesive or mutually attractive endothelial cells suffice for 
the formation of blood vessel sprouts and vascular networks. 
Because each mathematical representation of a model intro- 
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duces potential artifacts, it is important that model results 
are reproduced using alternative modelling paradigms. Here, 
we present a lattice-free, particle-based simulation of the 
cell elongation model of vasculogenesis. The new, particle- 
based simulations confirm the results obtained from the pre¬ 
vious Cellular Potts simulations. Furthermore, our current 
findings suggest that the emergence of order is possible with 
the application of a high enough attractive force or, alter¬ 
natively, a longer attraction radius. The methodology will 
be applicable to a range of problems in morphogenesis and 
noisy particle aggregation in which cell shape is a key deter¬ 
mining factor. 

Keywords Vasculogenesis • Cell-based model • Cell 
elongation • Morphogenesis • Alignment order 

1 Introduction 

The vascular system is one of the most important organs 
in large multicellular organisms. During embryogenesis the 
need for an efficient transport of nutrients and waste prod¬ 
ucts arises naturally as the increasing size of the organism 
makes diffusion less and less efficient. A similar situation 
arises in solid tumours, but instead of the de novo formation 
of the vasculature (vasculogenesis), tumours are able to re¬ 
model the vascular network in the surrounding tissues and 
promote formation of new branches from the existing ones 
(angiogenesis) [7]. Understanding the basic principles be¬ 
hind these processes could help in controlling tumour growth, 
as well as improving restoration of normal vasculature after 
trauma such as a stroke. It is tempting to speculate that the 
same basic principle that generally promotes network for¬ 
mation during vasculogenesis also drives blood vessel sprout¬ 
ing during angiogenesis. 

To explore mechanisms behind vasculogenesis it is ben¬ 
eficial to turn to computational modelling where all param- 
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eters are under control, before attempting experimental val¬ 
idation where unknown factors could complicate the inter¬ 
pretation of the results. Various hypotheses have been con¬ 
structed with the aid of computational biology [3], however, 
we argue that the best modelling approach to study network 
formation is through cell-based models [13]. Such models 
are well suited to describe collective behaviours emerging 
from cellular properties and interactions, such as collective 
cell motion [11], or various morphogenetic processes [14]. 
Vascular networks contain structures typically on the scale 
of a couple to tens of cells which would invalidate continuum- 
based approaches, that are better fitted for describing cellu¬ 
lar structures of at least hundreds of cells in size. 

Previous studies using cell-based modelling have identi¬ 
fied a handful of candidate mechanisms for sprouting and 
network formation [4]. By the nature of biological cells, 
all of these mechanisms include a form of long-range at¬ 
traction that ensures cell cohesion, and a short-range repul¬ 
sion that is needed to represent cell incompressibility and 
thus prevent unrealistically high cell densities. One class 
of models explains the long-range attraction via mechano- 
taxis, whereby cells contract and deform an underlying soft 
substrate, pulling cells towards each other [10] or where 
the substrate locally stiffens in response to cellular contrac¬ 
tion directing cells towards one another via durotaxis [16]. 
This mechanism suffices to explain how in some culture sys¬ 
tems, for example bovine aortic endothelial cells on poly¬ 
acrylamide gels [16], networks form only on substrates of 
specific stiffness. Contractile cells on deformable substrates 
have been described as force dipoles revealing that such dipoles 
tend to form locally or globally ordered structures depend¬ 
ing on the Poisson ratio of the substrate [1,2]. An alterna¬ 
tive mechanism, based on the observation that cells pref¬ 
erentially adhere to elongated cellular structures, explains 
how networks can form on bare substrates [19,20]. Another 
body of work focuses on chemotaxis towards a self-secreted 
chemokine, in combination with either contact-inhibition of 
chemotaxis [15] or cell elongation [12], that proves to be a 
robust mechanism albeit the chemokine remains to be iden¬ 
tified. While the different mechanisms of attraction distin¬ 
guish between different hypotheses, the mechanism of re¬ 
pulsion is less important and is either introduced explicitly 
(for example in partial-differential equation models [10] or 
particle-based models [20]) or is an implicit property of the 
model (as in the CPM [12,16] for instance). 

The model mechanism we focus on in this paper demon¬ 
strates that cell elongation and volume exclusion together 
with contact-dependent cell adhesion [17] or longer range 
cell-cell attraction [12], are sufficient for network forma¬ 
tion (Fig. 1). Whereas slightly adherent, dispersed cells ag¬ 
gregate into compact clusters (Fig. la), elongated cells form 
network-like structures in the models with (Fig. lb) or with¬ 
out (Fig. Ic) additional longer range chemotaxis. In these 



(a) (b) (c) 


Fig. 1 Cellular Potts simulations of chemotaxis and cell elongation 
models of vasculogenesis. a, A simple mechanism of chemotaxis to¬ 
wards a self-secreted chemokine generates cell aggregates, b, Chemo¬ 
taxis together with cell elongation gives rise to robust networks, c, Cell 
elongation without chemotaxis is sufficient to explain network forma¬ 
tion. Used with permission from [17] 


models of elongated cells, the branches are stabilised by the 
increasing rotational inertia of growing clusters of elongat¬ 
ing cells: a tightly packed branch of elongated cells is harder 
to displace or reorganise, therefore cells are frozen into the 
branched pattern in the model. A somewhat similar cellular 
process to elongation described by this model is the apical 
constriction of epithelial cells that plays a major role in, for 
example, gastrulation and neurulation [18]. Although api¬ 
cal constriction leads to cell elongation, this elongation may 
only be a passive result of volume conservation and cyto¬ 
plasmic fiow [8]. 

As each model implementation is a simplification with 
different limitations, it is worth to examine the model at 
hand using more than one implementation with different im¬ 
plicit assumptions associated with the simulation methodol¬ 
ogy. If the hypothesis is independent of the implementation, 
it should be able to reproduce the phenomenon in differ¬ 
ent implementations. For example, this has been performed 
in the case of the preferential adhesion hypothesis, which 
has been tested in both a lattice-free [20] and a grid-based 
model [19]. 

The previous studies [12,17] used the same cellular Potts 
model (CPM) description which could introduce model-specific 
artefacts. For example, highly elongated cells could breach 
the limitations of cell representation within the CPM, as cell 
elongation is achieved by maximizing the largest moments 
of inertia of the cell [12], which leads to unrealistic thicken¬ 
ing of the cell body at its extremities. As artefacts may result 
from the grid-based nature of the model and the implementa¬ 
tion of cell elongation, we chose to implement the cell elon¬ 
gation hypothesis in a lattice-free model using a generalized 
attraction and repulsion between cells, describing a class of 
models. We use this model to test if the hypothesis holds and 
to compare our results to previous reports. 
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Fig. 2 Cells and their interaction in the model, (a) Cells are repre¬ 
sented as 2D ellipses consisting of a repulsive core as cell body (blue). 
Cell attraction occurs at a finite range around the core in an area repre¬ 
sented as a bigger ellipse surrounding concentrically the repulsive core 
(red), (b) Cell-cell interaction is governed by the overlap area of the 
two inner ellipses (blue area) and the overlap of the two outer ellipses 
(red area) 


2 Methods 


2.1 Computational Model 


Each cell i in the model is described as an ellipse on a 2D 
plane, characterized by its position x/, direction of elonga¬ 
tion 0/, area A/, and the aspect ratio of major-to-minor axis 
Si. This ellipse is a repulsive core representing the incom¬ 
pressible cell body (Fig. 2a, blue) and is surrounded by a 
larger concentric ellipse that is responsible for cell-cell at¬ 
traction in the model (for example via filopodial adhesion) 
(Fig. 2a, red). Cell motion is modelled as a persistent dif¬ 
fusion process as in a previous study [20]. The change in 
velocity \i{t) for cell i at time t is described as: 


dyi{t) _ J_ 

dt Mi 







( 1 ) 


Here mt is the mass of cell i (mi = m = 1), a is a damping 
parameter, is a. uniformly distributed random vector in 
the 2D plane e [—0.5 : 0.5] x [—0.5 : 0.5]), and model 
parameter Ny sets the amplitude of this translational noise. 
The last term is the sum of pair interactions between cells 
describing short-range core repulsion, and a long-range but 
finite attraction: 


Fij = XrAr - Xa^a^ ( 2 ) 

Ar is the overlap area between the two repulsive ellipses of 
the two cells (Fig. 2b, blue area), and Aa is the overlap area 
of the two outer ellipses (Fig. 2b, red area), and Xr and Xa 
are model parameters. The overlap is calculated using a pre¬ 
viously published method [9]. 

Model cells are rotated to minimize overcrowding and 
maximize attraction [2]. For every interval Ar an attempt 


is made to change the orientation of N randomly selected 
cells with a random angle where G [— TT : ;r] is 

a uniform random variable. The change is accepted with a 
turn-probability 


: mm 


L I ^ . 


( 3 ) 


where F^j is the pair interaction in a configuration after the 
proposed rotation and Na is the angular noise parameter. 

Simulations were initiated with N = 1000 overlapping 
cells distributed at random within a 200 fimx 200 fim unit 
area. Cell area was fixed at 250 /im^ throughout the study. 
The size of the outer “adhesive” ellipse is set by the at¬ 
traction radius parameter defined as Ra = aaja, where Ua 
is the semi-major axis of the adhesive (outer) ellipse and a 
is the semi-major axis of the repulsive (inner) ellipse, and 
the aspect ratio of the outer and inner ellipses is the same. 
Baseline parameters used are: Ra = 1.5, Xa = 0.0006, = 

0.02, a = 0.4,A^v = 0.4, Na = I and s = 7. These parame¬ 
ters were estimated by macroscopic inspection of the model 
to yield “realistic” pattern formation as a measure of val¬ 
idation. To prevent excessive cell overlap, the magnitude 
of the repulsion strength (Xr) is two orders of magnitude 
larger than the magnitude of the attraction strength (Xa). The 
model’s sensitivity to the parameters was tested by alter¬ 
ing the parameter values until unrealistic results, either con¬ 
cerning cell movement or cell interactions, were produced. 
Model equations 1 and 2 are integrated numerically using 
the forward Euler method with a fixed time-step of At = 0.1 
(or At = I for the longer simulations of up to r = 100,000 
shown in Fig. 5), and orientations are updated after each it¬ 
eration synchronously using Eq. 3. Integration was stopped 
at r = 25,000 (or ait = 100,000 for simulations shown on 
Fig. 5). In each case the system reached a quasi-stationary 
state at the end of the simulations. An implementation of this 
particle-based system is provided in Online Resource 3. 


2.2 Order parameter 

Focal alignment of elongated cells has been shown to play 
an important role in the formation of networks previously 
[17]. To measure alignment, a local orientational order pa¬ 
rameter is used: 

S(r) = {cos(2ei(r)))i. (4) 

6i(r) is the polar angle between (j)i, the orientation of cell 
/, and rii(r) = ( 07 )|y |x -x |<r}’ average cell orientation 
within r distance of cell i. S(r) is isotopic and takes the value 
of 0 for randomly oriented cells and 1 for perfectly aligned 
cells. Smaller values of the order radius r describe the align¬ 
ment in the close vicinity of the cell (r = 20), or in the local 
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Fig. 3 Ranges of local order. Two values of r for the order parameter 
calculation shown for the centre cell within a typical cell arrangement. 
While r = 20 describes order at the local cell neighbourhood, r = 80 is 
an indication of order in multicellular structures in the vicinity 

structures (r = 80, Fig. 3), while any global order is captured 
by the global order parameter, Sg = S{r oo'), where all the 
cells are considered for the calculation of r]/. The software 
used for measuring the order parameter is provided in On¬ 
line Resource 3. 


3 Results 

3.1 Elongated cell shape induces network formation 

Rounded cells (s = 1.5) aggregate into less ordered clus¬ 
ters (Fig. 4a and the animation in Online Resource 1), while 
more elongated cells with aspect ratio s = 1 in the simula¬ 
tion interconnect to form elongated branches and networks 
(Fig. 4b and the animation in Online Resource 2). The time 
evolution of the order parameters calculated for two simula¬ 
tions of round and elongated cells shown in Fig. 4 demon¬ 
strates that the initially disordered cells align within t=10^ 
and keep a similar level of order for at least 10 times longer 
at finite length scales (5(20) and 5(80) on Fig. 5a). Positive 
value of the global order parameter (Sg) in the elongated cell 
configurations results from finite system size and indicates 
the emergence of a system-wide alignment (Fig. 4c). Elon¬ 
gated cells (s = 7) produce a significantly higher order than 
more round cells (s = 1.5) at all length-scales (Fig. 5b). Or¬ 
der at the cell-to-cell range (5(20)) always supersedes order 
at the multi-cellular level (5(80)), which is always higher 
than global order (Sg). 

Simulations with cells of increasing aspect ratios and 
constant areas reveal that order is gradually increased in sim¬ 
ulations with more elongated cells (Fig. 6). When simula¬ 
tions with increasing 5' are compared, order emerges first at 
short length scale (5(20)) followed by longer range order. 



(c) 


Fig. 4 Cell elongation induces network formation a, Short cells with 
aspect ratio .y = 1.5 form compact aggregates, b, Elongated cells with 
aspect ratio of .y = 7 form networks, c, Temporal evolution of network 
formation for cells with aspect ratio of .y = 7. Global order is estab¬ 
lished at about t = 10^. The system reaches a fluctuating steady state: 
the morphological details change over time, but the morphological fea¬ 
tures remain statistically constant. 

Global order appears in simulations with s > 4. The high 
variation of the global order among the different simulations 
is consistent with the fluctuation shown on Fig. 5a. 


3.2 Increased range of cell attraction enhances cell 
alignment 

Previous studies showed that chemotaxis could play an im¬ 
portant role in vascular network formation 112, 15]. Since 
chemotaxis may be interpreted in our model as an adhe¬ 
sive interaction over a longer range, we investigated how 
the range of attraction affects the observed cell alignment. 

As expected, a short range attraction (Ra = 1.1) results in 
a dissociated cell configuration and low order (Fig. 7a and 
Fig. 7b, Ra = l.l), compared to the control parameters (Fig. 4b). 
A long attraction range, however, results in a marked in¬ 
crease in ordering, consistent with the longer range attrac¬ 
tion of the chemotaxis studies. Interestingly, global order in¬ 
creases markedly at an interaction range of less than twice 
the cell size, where long, aligned strands of cells are formed 
(Fig. 7d, Ra = 2). 
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Fig. 5 Order parameters of elongated and rounded cell simula¬ 
tions. a, Temporal evolution of the order parameter, S{r), for r = 20 
(green curves), r = 80 (red curves) and r^oo (blue curves) for two sys¬ 
tems of more rounded cells (s = 1.5, dotted lines) and elongated cells 
(s = 7, solid lines; At = 1). Ribbons show standard deviations over 10 
runs. Finite scale order reaches a plateau at t = 10^ and remains ap¬ 
proximately stable up to at least t = 10^. b. Order parameters from 10 
independent simulations show that alignment at ^ = 10^ is higher for 
elongated cells than rounded cells at all length scales {At = 1). Within 
configurations order at shorter scales is always higher than order at 
longer scales 


Increased alignment is also achieved by increasing at¬ 
traction strength through parameter (Fig. 8), while in¬ 
creased repulsion leads to dissociation and consequent loss 
of order (Fig. 9). 


3.3 Noise is not essential for alignment 

We assessed the importance of noise for cell alignment in 
our model (Fig. 10). In the absence of translational noise 
(Fig. 10b) or low rotational noise (Fig. lOe) cells align and 
form networks. At high noise levels the two types of noise 
act in similar manner; high translational noise results in a 
decay of global order (Figs. 10a and 10c), while high angu¬ 
lar noise allows for more energy inefficient rotations, thus 


1 . 00 - 



Fig. 6 Emergent order in the model is an increasing function of cell 
elongation, a. Simulations with uniform cells of various aspect ratios 
show that slightly elongated cells align locally, while more elongated 
cells give rise to order at longer length scales. Ribbons show standard 
deviations over 10 runs after t = 2.5 x 10^. b-d. Configurations ait = 
2.5 X 10^ for aspect ratios .y = 1.5, .y = 4.5, and s = 7 



Fig. 7 Increasing range of attraction increases order, a, Local, in¬ 
termediate, and global order are higher with larger attraction ranges 
for cells with aspect ratio of 5 = 7. Ribbons show standard deviations 
over 10 runs after f = 2.5 x 10^. Configurations at f = 2.5 x 10^ for (b) 
short, (c) medium and (d) long attraction ranges Ra = 1.1, = 1.6 

and Ra = 2 
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Fig. 8 Stronger attraction leads to increased order, a, Order param¬ 
eters from simulations with increasing attraction parameters for cells 
with aspect ratio of ^ = 7. b, Simulation configurations with low at¬ 
traction strength = 4 x 10“^. c Simulation configurations with high 
attraction strength = 2 x 10“^ 



(a) 












(b) 

(C) 


Fig. 9 Repulsion acts against order, a, Order parameters from simu¬ 
lations with increasing repulsion parameters for cells with aspect ratio 
of ^ = 7. b, Simulation configurations with a low repulsion strength of 
Xr = 10“^; c, Simulation configurations with a high repulsion strength 
of A, = 9 X 10-2 


1 . 00 - 



o!o o!5 i!o i!5 2.0 

Angular Noise Factor 


(d) angular 



(e) (f) 

Fig. 10 Effect of noise on alignment, a, Order parameters as a func¬ 
tion of vectorial noise for cells with aspect ratio of 5 = 7; b, Simulation 
configurations in absence of translational noise; c, Simulation configu¬ 
rations with increased translational noise {Ny = 3); d, Order parameters 
as a function of angular noise for cells with aspect ratio of ^ = 7; e, Sim¬ 
ulation configurations with low angular noise {Na = 0.1); f, Simulation 
configurations with increased angular noise (Na = 2). Order parameters 
and configurations are shown at f = 2.5 x 10^ simulation time 
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(a) (b) (c) 

Fig. 11 Effect of parameters on alignment without cell elongation. 

a, Simulation configurations with increased radius of attraction Ra=2\ 

b, Simulation configurations with increased attraction strength Xa = 

2 X 10“^; c, Simulation configurations with low repulsion strength Xr = 
10“^. Configurations shown at f = 2.5 x 10^. All cells have an aspect 
ratio of .y = 1.5 

decreasing global order (Figs. lOd and lOf). At low noise 
levels, the two types of noise act differently; while low trans¬ 
lational noise has no or very little effect on the orientational 
order (Fig. 10a), low angular noise results in higher local or¬ 
der and lower global order (Fig. lOe), suggesting the forma¬ 
tion of isotropic branches that point in different directions. 

3.4 Cell elongation is essential for alignment 

Finally, we tested whether cell elongation is required for 
alignment, by influencing the parameters that had positive 
impact on cell order in the previous simulations. Using rounded 
cells with aspect ratio of 5' = 1.5, three scenarios with in¬ 
creased range of attraction (Ra), increased attraction strength 
(Xa) and reduced repulsion strength (Xr) were tested. Or¬ 
der parameters from these simulations together with order 
parameters from simulations of round and elongated cells 
(Fig. 4a-b) are summarized in Table 1. These indicate that 
neither increased adhesion strength or range, nor decreased 
repulsion is able to order the cells to an extent compara¬ 
ble with elongated cell simulations. Without elongation, cell 
clusters are unable to deviate from their compact aggregates 
(Fig. 11). 


Table 1 Order parameters for round (s = 1.5) and elongated (s = 1) 
cells, compared with cells with aspect ratio of 5 = 1.5 at f = 2.5 x lO'^ 
tested with parameters that favour cell ordering 


Simulation 

Parameter 

5(20) 

5(80) 

5'(inf) 

round cells 

5= 1.5 

0.38 

0.13 

0.03 

elongated cells 

S = 1 

0.88 

0.67 

0.23 

high attraction range 

(N 

II 

0.31 

0.11 

0.04 

high attraction 

Aa = 2 X 10-3 

0.38 

0.08 

0.02 

low repulsion 

II 

o 

0.36 

0.12 

0.03 


4 Discussion 

A simple particle-based model is introduced to demonstrate 
that elongation can indeed aid cells to aggregate into a net¬ 
work. The mechanism for network formation is simple: Cell 


elongation, together with attraction and core repulsion, leads 
to local cell alignment and increase in orientational order. 
The fact that this prediction holds in two completely differ¬ 
ent model implementations, our particle-based model and 
the CPM [12,17], gives confldence that the predictions are 
due to the explicit model assumptions, not due to unintended 
properties of the simulation methodology. Time evolution 
of order and dependence on elongation (Fig. 5b) are simi¬ 
lar in both implementations. Increasing the attraction range 
of cells in the current implementation allows the study of 
the transition between the pure elongation hypothesis [17] 
to the chemotaxis and elongation hypothesis [12]. Extended 
attraction range, in the form of chemotaxis, has been shown 
to help in the formation of regular networks with a more 
deflned pattern size across the whole system [12]. Consis¬ 
tently with this previous result, we found that a longer range 
of attraction results in higher order and a more network-like 
structure (Fig. 7). 

Cell movement dynamics is a marked difference between 
the CPM and our model. In the CPM cell movement emerges 
from the displacement of the cell boundaries, yielding a more 
amoeboid cell movement. By contrast, in particle-based mod¬ 
els model cells translate as a single unit, including translo¬ 
cation and rotation in our case. Cell movement in the model 
is described by an overdamped dynamics, where the damp¬ 
ing factor (a in Eq. 1) controls the movement persistence of 
cells. For large values of a, the left-hand-side of Eq. 1 van¬ 
ishes and Eq. 1 can be approximated by an algebraic equa¬ 
tion giving the cell velocities as a function of the external 
forces acting on the cells. We have chosen here to work with 
a differential equation description, thus keeping the damp¬ 
ing factor in our model as an explicit parameter. Although 
we have here only studied overdamped kinetics, smaller val¬ 
ues of a could mimic persistent cell motility: it takes time 
for a cell to change its direction if, e.g, the chemoattractant 
gradients change; the actin cytoskeleton needs to reorganize 
which takes time. Such persistent motion was also explicitly 
included in early chemotaxis-based partial-differential equa¬ 
tion models of vascular network formation [5]. Although 
persistence of motility was later shown to be unnecessary for 
network formation, it was argued that in the CPM the cell 
shape may cause some directional persistence (Figure 11 
of Ref. [12]). In the present work directional persistence 
is severely reduced and is largely independent of the cell 
shape, suggesting that persistence is not required for net¬ 
work formation. 

Stochasticity in our model is introduced through the an¬ 
gular noise similar to the noise of the Vicsek model [21] and 
the translational noise of the Gregoire - Chate model [6]. 
Note that these noise factors describe fluctuations in the sys¬ 
tem, and therefore the considered noise factors are inde¬ 
pendent from one another and from the noise in the past. 
Our results from the study of the angular noise show that 
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noise hinders the global order, in good agreement with pre¬ 
vious particle-based simulations using only angular noise 
and force-dipoles to describe multicellular structure forma¬ 
tion [2]. 

Interestingly, both an increase in attraction range (Fig. 7) 
and attraction strength (Fig. 8) leads to a higher tendency to 
global ordering as packing of the cells becomes tighter. This 
is similar to the previous observation in the CPM models 
[12,17], where the more adhesive cells become more com¬ 
pact and highly ordered in domains, but without the emer¬ 
gence of global order. In the CPM cell shape is an emergent 
property of the cells and therefore cells are able to deform 
in order to maximise the contact surface and compactness 
(see Fig. 4a of [17]). In contrast, cells in our implementation 
are unable to deform such that they might be less able to ac¬ 
commodate the boundaries of such locally ordered domains. 
This may force them to align globally through the compact¬ 
ing force of the strong attraction. Secondly, when attraction 
range is increased to a distance of two cell diameters in the 
current implementation, alignment order spans through do¬ 
mains of at least 1000 cells, with strands of cells aligned 
in parallel even without tight packing (Fig. 7d, 7?^ = 2). An 
effect contributing to this global ordering may be the im¬ 
plementation of noise in the current, particle-based imple¬ 
mentation. In the CPM, at high densities, cells hinder each 
other’s translation and rotation, resulting in an increasingly 
slow development of the pattern as the branches grow [17]. 
By contrast, in the particle-based model translation is not 
hindered by adjacent cells and occurs for cells as a whole 
(Eq. 1), while cells continue to rotate independently as a 
whole as long as conflicts with adjacent cells persist (Eq. 3). 

In conclusion, here we introduced a particle-based model 
to re-examine a hypothetical mechanism for the formation 
of micro vascular networks: i.e., that elongated vascular cells 
tend to aggregate into branches of a network structure. Our 
previous work [12,17] simulated this potential mechanism 
using the cellular Potts model, demonstrating that indeed 
elongated cell shape, in combination with mutual attraction 
or adhesion suffices for the formation of network-like pat¬ 
terns. Here we have shown that this phenomenon also occurs 
in a lattice-free particle-based model, adding confidence that 
the effect is not caused by artifacts of the cellular Potts model 
or of the present, particle-based model. Thus our models 
suggest that network formation is a natural emergent prop¬ 
erty of elongated, adhesive objects in a stochastic system. 
Because a range of alternative mechanisms for network for¬ 
mation and angiogenic sprouting have been suggested (re¬ 
viewed in [14] and [3], see also [16]), to what extent the 
present mechanism contributes to angiogenesis in vivo or 
in vitro at this point must remain an open question. These 
and similar studies, however, help to generalize and catego¬ 
rize the main requirements for network formation and an¬ 
giogenic sprouting, thus contributing to ongoing efforts to 


identify the controlling factors of angiogenesis from a bio¬ 
physical point of view. 

Online Resource 1 Simulation video of rounded cells for 
r = 2.5 X 10^. Simulated cells of aspect ratio s = 1.5 form 
aggregates. Every frame is after t = 50. Every second of 
video is ^ = 1250; see https: //youtu. be/BgQwdlaGxAI 
Online Resource 2 Simulation video of elongated cells for 
t = 2.5 X 10"^. Simulated cells of aspect ratio s = 1 form 
networks. Every frame is after r = 50. Every second of video 
ist = 1250; see https://youtu.be/9zUGEBk6pio 
Online Resource 3 implementation of the presented 
particle-based system, and software for measuring the order 
parameter 
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